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ABSTRACT 

An antenna shower is a parton shower in which the basic move is a 
color- coherent 2 — )■ 3 parton sphtting process. In this paper, we give com- 
pact forms for the spin-dependent antenna sphtting functions involving 
massive partons of spin and spin ^. 
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1 Introduction 



The modeling of physics at high energy colhders rehes heavily on our under- 
standing of QCD. Quarks and gluons - collectively, partons - that are produced in 
high energy reactions are observed as jets of hadrons. The structure of each jet is 
determined by the pattern of radiation of additional partons from the original one 
produced in the central hard scattering reaction. For this reason, much attention 
has been given the past few years to the development of methods for creating parton 
showers, systems of partons created with the distributions predicted by QCD. 

The traditional approach to the generation of parton showers is based on splitting 
off partons through a 1 — ?■ 2 branching process. This philosophy is incorporated 
in the widely used event generator programs PYTHIA [1] and HERWIG [2]. The 
construct of building a shower from 1 — )■ 2 branching, often called a 'dipole shower', 
omits an important aspect of the physics. The longitudinal momentum distribution 
in the 1 — ?■ 2 splitting is given by the Altarelli-Parisi splitting functions In QCD, 
partons are emitted coherently from the two legs of a color dipole. The emission 
amplitude is then enhanced inside the dipole and, more importantly, cancels outside 
the dipole. In the 1980's, Marchesini and Webber argued that this effect could be 
incorporated into dipole showers by imposing angular ordering of emissions Thus, 
HERWIG is built around an angular-ordered parton shower, and PYTHIA, though 
it uses a different ordering scheme to choose its branchings, vetos emissions that are 
out of angular ordering. 

Alternatively, one might build up a parton shower directly from the color dipoles, 
using the 2 — )■ 3 process of emission of a parton by a dipole as the basic branching 
process. This construct is called an 'antenna shower'. The scheme was realized in the 
program ARIADNE, by Andersson, Gustafson, Lonnblad, and Pettersson ^ and, 
more recently, by the program VINCIA, by Giele, Kosower, and Skands [6j. The 
approach is of interest both in creating new parton shower codes for the purpose of 
matrix element-parton shower matching and because of its promise to yield a more 
accurate treatment of color dynamics in parton showers. 

Recently, there has been much interest in the tagging of boosted heavy particles 
such as the top and Higgs observed as exotic jets [7]. Since tagging methods rely 
heavily on color flow, it is interesting to have a variety of approaches to the simulation 
of color flow in parton showers in order to test the robustness of these algorithms. 

We have been engaged in providing a well-defined foundation for antenna show- 
ers, giving explicit calculations of the splitting functions that generate these showers 
and generalizing previous work to spin-dependent formulae. In a previous paper, 
we presented the complete set of spin-dependent antenna splitting functions needed 
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to describe quark and gluon parton showers In this paper, we continue our 
study of this approach by presenting the spin-dependent antenna sphtting functions 
for showers with massive particles. In constructing a shower for massless particles, 
spin-dependence is a convenience, especially for matching with full QCD amplitudes. 
For massive particles, it is more important to preserve spin information, because the 
decays of heavy particles such as the top quark are spin-dependent and so the exper- 
imental acceptance for the heavy particles varies significantly with their longitudinal 
polarization. 

The formalism presented here has the same strengths and weaknesses as our pre- 
vious work. We will calculate in the kinematics of final-state showers, using effective 
operators of definite spin to represent the 2-particle color dipole state before the split- 
ting. We will work in the limit of a large number of colors in QCD for which the 
concept of a color dipole is strictly defined. Within this approximation, we will derive 
formulae for splitting functions with any ratio m/Q between the mass of the particle 
and the mass of the two-particle system. These formulae will necessarily be less sim- 
ple than those found in [8] for the massless case. We will see, though, that we can 
make use of spinor product formalism |9j to write these splitting functions relatively 
compactly. The simplicity of these expressions is connected to their relation to the 
Maximally Helicity Violating amplitudes of QCD. This point was originally made for 
the massless case in [10] and is discussed in some detail in [8j. 

The formalism of QCD antennae was originally developed as a tool for the sub- 
traction of infrared divergences in higher-order QCD calculations. This approach 
to QCD calculation was pioneered by by Kosower [TT|T2] . Gehrmann-De Ridder, 
Gehrmann, Glover, and their students have developed this approach into a sophisti- 
cated method applicable to NLO and even NNLO computations [T^P^ . Using this 
formalism, Gehrmann-De Ridder, Gehrmann, and Glover have proposed forms for the 
spin-summed antenna splitting functions of massless quarks and gluons [T5|IT6] . Our 
previous paper reviews this latter work and compares the results from our method to 
theirs. There is no universal form for antenna splitting functions. The behavior of the 
splitting functions is prescribed in the soft and coUinear limits but, away from those 
limits, different expressions are possible, depending on the framework used in the 
derivation. The systematic differences between the different proposals are explored 
in [8]. 

Following the methods of P^IITU] . splitting functions for massive, spin summed 
antennae were constructed in [T7|18fl 9]. Again, our expressions agree with these 
in having the correct soft and quasi-collinear behavior but differ away from these 
limits. The addition of mass greatly complicates both the expressions for the splitting 
functions and the precise specification of the boundaries of phase space. Because of 
this, we do not present a detailed comparison to other massive splitting functions 
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here. 

The outhne of this paper is as follows: In Section 2, we will analyze the case 
of gluon radiation from an antenna composed of a massive spin-| fermion (Q) and 
a massless spin ^ fermion (g) in a configuration of zero helicity. All of the new 
complications that arise when we deal with massive particles can be illustrated in 
this context. We will write expressions for the splitting functions in terms of spinor 
products of lightlike vectors associated with the massive vectors of the particles before 
and after the splitting. In Section 3, we will discuss the kinematics of these massive 
splittings and the evaluation of the the spinor product expressions. 

With this introduction, we can go systematically through the various cases of 
antennae composed of massive and massless particles. In Sections 4 and 5 we will 
analyze in turn the cases of antennae with spin and spin | massive particles recoiling 
against quarks and gluons in which the antennae emits another quark or gluon. In 
Section 6, we discuss the analysis of the general case of a pair of massive particles, 
spin or spin |, radiating gluons. In Section 7, we discuss antennae that create a 
pair of massive particles. Section 8 gives some conclusions. We collect the complete 
set of massive antenna splitting functions derived in this paper in Appendix A. 



2 The spin zero fermion-quark antenna 

The simplest case of a splitting function with massive particles arises in the system 
of a massive and a massless fermion created by a spin operator. In this section, 
we will work out the spin- dependent splitting functions for this case following the 
prescriptions in [8]. We will then discuss the interpretation of these formulae and 
their comparison to the standard Altarelli-Parisi splitting functions for a massive 
quark [2D]. 

In [S] , each case of a spin-dependent splitting is associated with a gauge-invariant 
operator that creates the antenna. For this case, the required operator is 

O = QlQr (1) 

where q is an ordinary quark whose mass can be ignored and Q is a massive quark. 
This operator creates a 2-particle state 

QlQl (2) 

with total spin zero about the production axis. Antennae with overall opposite helicity 
or with antiquarks have the same splitting functions, by the P and C invariance of 
QCD. 
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In [8], we wrote the basic formula for final-state antennae splitting of massless 
particles in the following way: Notate the splitting as AB — )■ acb, with 

{A + By = SAB = Q' . (3) 

Throughout this paper, for any 4- vectors i, j, we will define 

Sij = {i+jY = m^,+2i-j+m^j. (4) 

Let Za, Zh, Zc be the momentum fractions of a, b, and c relative to their maximum 
value, 

Za = —^jT' ' Za + Zb + Zc = 2 . (5) 

Then the probability of a splitting is given by 

J ^P^ot) = A^c^(^) J dZadZbS{Za,Zb,Zc) • (6) 

where A^c = 3 is the number of colors in QCD and K is the momentum of the partons 
in the center of mass system of the original 2-particle antenna. In the massless case, 
Q/2K = 1. The distribution S is the splitting function. In |8], we computed this 
function as the ratio of 3- to 2- body amplitudes of an appropriate local operator. 



S = Q' 



M{0^acb) 



MiO^ AB) 



(7) 



This formula is still correct for the massive particle antennae discussed in this paper. 
We will discuss the kinematics of these antennae in more detail in Section 3. 

In the limit in which c becomes coUinear with a or b, the antenna splitting functions 
reduce to the Altarelli-Parisi functions P{z) that describe 1 — 2 splittings. For this 
limit, the formulae are not as simple in the massive case as they are in the all-massless 
case. We will present the explicit formulae and check them for the spin zero antenna 
later in this section. 

To compute the amplitudes in ([T]), we use the spinor product formalism for massive 
particles of Schwinn and Weinzierl [21]. For a massless particle, the states of definite 
helicity are well-defined and Lorentz invariant. For a massive particle, the spin states 
depend on the frame chosen to evaluate them. In the Schwinn- Weinzierl formalism, a 
massless reference vector q is used to define that frame. The spinors for an outgoing 
massive fermion of mass m are written 

_ [q{p + m) _ {q{p + m) 

MP) = uM = , (8) 



4 



where the flatted vector is defined by 



2 

m 



p =p - q (9) 

2q ■ p 

A particularly useful choice for q is the lightlike vector in the opposite direction from 
p. Rotating coordinates so that 

p={E,0,0,p) with = + , (10) 

let 

p« = ^(E + p)(l, 0,0,-1) , (11) 
Then if we set q = p^, the fiatted vector is 

p' = ^{E+p){l,0,0,l) . (12) 



This is very convenient. With this choice of q, the spinors defined in ([8j) are just 
the usual spinors of definite helicity. Using the basis of Dirac matrices where 7^ is 
diagonal, it is easy to see that (|8|) reduces to 



^L = (v^ ^)®(o 1) 

^R={^/W 0) • (13) 



Using these conventions, we can easily compute the 2 particle matrix elements 
of the operator ([T]). Denote the momenta of the initial-state heavy quark and light 
antiquark as A and B, respectively. Then 

MiQm.) - ^ (14) 

The helicity of the q must be L, but the heavy quark created by ([T| could be in 
either spin state. However, with the usual definition of helicity, the production of 
Qr(Il from a spin operator would be forbidden by angular momentum. Indeed, 
when we set q = A^, 

M{Qru) ~ (A^B) = , (15) 

because Ai is a lightlike vector parallel to B. The only nonzero matrix element is 
then 

MiQLU) = {A'B) ; (16) 
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this gives the denominator in (|6|. It is convenient that 

\{A'B)\' = Q''-m^ = 2QK , (17) 

with i^' as in ([6]). 

It is straightforward to work out the numerator of (|6| for the four possible spin 
states of the 3-particle system Qgqi- As in [8J, we label the three final-state momenta 
as (a, c, b), with the emitted particle as c. The results, using a general reference vector 
g in ([8|, are 



MiQugLQi) 
MiQugnqL) 



1 f {ca%Qb) ^ [qQa') 
Sac - rn? [be] 
a!'b)[c( 



[qc] I 



{bc){sac - m?) 



m 



f {cq)[a'Qb) ^ [a'Qq) \ 



[a^c]{qa^) I Sac - m 

m{qb)[cQb) 
{bc){qa^){sac - m? 



[be] J 



We have omitted the overall factor of {gT"-)- When we put q = a*, we can recognize 
the simplification 

[a^Qa^) = [a^Qa^) = . (19) 

This follows from the fact that the 4-vector Q is a linear combination of the two 
lightlike vectors and a". Now square these expressions and combine with (17) to 
evaluate ([7|. This gives 



SiQigLqi) 

SiQEgmLi 
SiQugLqi) 
S{QRgRqL) 



Q_ 
2K 
Q_ 
2K 



(a^c)[at 



[af c] [cac) 
{a^b)[cQb) 



{be) [cac) 





{a^c)[d'Qb) 


2K 


{a^a^)[a^c][cac) 


w?Q 


{a%)[cQb) 


2K 


{a^a^){bc) [cac) 



(20) 



In the all-massless case, we managed to produce antenna splitting functions that 
were simple rational functions of the Za pj. Here, the antenna splitting functions 
are more complicated, but not excessively so. The main complications come from the 
denominators {sac — fn^) = [cac), which do not factorize simply, and from the multiple 
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lightlike vectors needed to characterize the state of the massive quark. In this case, 
it is not so difficult to write the splitting functions in terms of 4-vector products: 



Qs^,,{2J-Qb-Q-a^-bQ^) 

QsM2b-Qc-Q-b-cQ^) 
K Sbc{sac - m?y 
Tn'^Qsj,{2a^ -Qb-Q-a^ -bQ^] 



c/^ -N m^Q SMi^b- Qc- Q -b- cQ^) 



However, the structure of the expressions is more clearly visible in the form (20). 



The expressions ( 20 ) contain exact tree-level matrix elements for the transition of 
the operator O to a three-particle state. They are correctly used in a parton shower 
for any values of m/Q and Pt/Q among the final-state particles, as long as the 
virtuality at the previous and successive branchings of the shower are well separated 
from Q. In the all-massless case discussed in [S], we made approximations to the 
splitting functions valid in the soft and coUinear limits. It is less obvious here which 
approximations are appropriate, and, in any case, we did not see how to achieve much 
further simplification. So we will stop at this point for this set of splitting fuctions 
and for all of the massive particle splitting functions quoted in this paper. 



To evaluate expressions of the type of (20), we find it easiest not to convert the 



expressions in (20) into 4-vector products or dimensionless scalars built from these 
but, rather, to directly evaluate the spinor brackets. We will discuss a strategy to 
evaluate these brackets in the next section. 

Finally, we must discuss the collinear limits and the connection to the the Altarelli- 
Parisi splitting functions. For the spin zero antennae, this connection is easiest to 
discuss for the limit c \\ b, where only massless particles are involved. We must still 
take account of the fact that, because b and c recoil against a massive particle, their 
maximum momentum is limited. To account for this, let 

'^"''^ " (1 - m7g2) ■ (22) 

so that Zb and Zc run from to 1 and, in the limit c \\ b, Zb + Zc = 1- Then, in this 
collider limit, S has the singularity 

S ~ 5a,A — PB^c{Zc) . (23) 



Sb, 



•c 
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The expressions in (20) satisfy this relation. The sphtting functions to QRgL,R niust 
have no coUinear singularity. This follows from the fact that [a"Qb) and {a^b) vanish 
when b becomes opposite to a. The cases of QLgL,R do have singularities proportional 
to s~,^ and s^^^, with the correct coefficients to match (22). 

In the limit c || a, where the 1 — ?■ 2 splitting involves a massive particle, the 
limit is slightly more complicated. For the splitting of a massive particle, the usual 
Altarelli-Parisi formula for the collinear splitting is conventionally rewritten as 

<<Prob = iV.|/d./|^;^M^P(..p.). (24) 

We divide the usual expressions for P{z,pt) by 2 so that these functions give the 
contribution from one of the two antennae that contribute to a collinear singularity. 
Mass-suppressed terms can contain an additional factor of (p^ + z^m?) in the de- 
nominator; this is why we have allowed the Altarelli-Parisi function to depend on px- 
With this formalism, for c becoming parallel to a, 

g2 

S{Za,Zb,Zc) ^P{zc,Pt) (25) 

Sac ^ "^yl 

where Sac = (a + c)^. Here again, the parameter Zc must be scaled to equal 1 at its 
maximum value, as in (22). For the present case in which the (ac) system recoils 
against a massless parton, Zc = Zc. 

To discuss the limits c || a, we first need to recall the Altarelli-Parisi functions for 
splitting of a gluon from a massive fermion. The Altarelli-Parisi functions are defined 
in the limit of not only collinear but also high energy emission. For a particle of 
energy E splitting to particles with transverse momentum pt and finite masses rrii, 
these functions describe the regime pr ~ "^i ^ E. For a splitting Q — )■ gQ, as we 
have in this case, the spin-summed splitting function is [20] 



P(z) = ^ 26 

z a ■ c 

This expression becomes clearer when it is written as a set of spin-dependent Altarelli- 
Parisi functions. In the convention defined by (24), 

p{Ql ^ Qigi) ^ 



p{Ql ^ QLgn) 



Pj^ + z'^m^ z 

V\ (1 - zf 
p\ + z'^m^ z 

1 4 

m z 



p{Ql ^ QngL) = , ^ , , 

Pj. + z^m^ z 

p{Ql ^ QRgR) = (27) 
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The sum of these terms does reproduce (26). The placement of the factors of z 
implements the dead cone in which soft radiation from a massive particle is suppressed 
within a cone of size I/7, where 7 is the boost of the heavy particle 



We can now compare the c || a limits of our antenna splitting functions to (27). 
In the collinear limit, 

...-m^-P^i^. (28) 

z[\ - z) 



Using this formula and the collinear limits of the spinor products, we find that (20) 
does satisfy (25) with (27), up to corrections of relative order m? /Q"^. In particular. 



in the limit c || a, becomes collinear with h. Then the vanishing of {a%) with no 



compensatory vanishing in the denominator gives the zero in the last line of (27). 



The spin-dependent splitting functions in the remaining sections of this paper also 
satisfy these checks on the collinear limits. For convenience, we list the complete set 
of mass-dependent, spin-dependent Altarelli-Parisi splitting functions that are needed 
for these checks in Appendix B. 



3 Kinematics of massive antennae 



The splitting functions computed in the previous section were written in terms 
of spinor products of massless vectors associated with the massive 4-vectors of the 
antenna. One should ask, how are these massless vectors computed? A similar 
question arises in the context of the formula (|6| for the antenna splitting probability. 
This equation is easily written down as the ratio of a cross section to produce a 3- 
body final state, integrated over 3-body phase space, to the cross section to produce 
a 2-body final state, without a radiated parton, integrated over 2-body phase space. 
In particular, the integral / dzadz^ is an integral over 3-body phase space. One should 
ask, what is the boundary of the region of integration for these variables, and how 
does one sample points in the interior of this region? 

For massless antenna, the answers to these questions are straightforward. For 
antenna with both radiators in the final state (FF antennae in the notation of [S]), 
the complete phase space region is the triangle 

<Za,Zb<l Za + Zb>l (29) 

and the region well described by the radiation process AB acb, with c soft, is the 
smaller region where 

< Zc < Za < I and < z^ < z^ < 1 (30) 
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To create an additional radiated particle in a state with massless particles, we 
choose a color-connected pair of particles AB, boost so that A and B are of equal 



length and back-to-back, choose {za, z^) as a random point in the region (30), replace 
the 2-particle system AB by the chosen 3-particle system acb, and, finally, reverse 
the boost to bring acb back into the original frame. The corresponding phase space 
regions and algorithms for antennae including initial-state particles are described in 
[8]. In this paper, however, we will only discuss final-state showers. 

We believe that these 4- vector configurations for massless particles provide a good 
starting point for constructing 4- vector configurations that include massive particles. 
Given a point {£i} in the phase space of of N massless particles, one can obtain a 
point {ki} in the phase space of N massive particles by rescaling 

ki = Xl (31) 

where A obeys 

^4 = i?CM, With 4 = (|A£l|2 + m2)i/2. (32) 

i 

Conversely, every point of the massive phase space can be constructed uniquely in 
this way. The scale factor A is close to unity unless one of the massive particles is 
nonrelativistic. The relation of the phase space measures for the massive and massless 
variables is |24] 

We will refer to the massless vectors as the backbone of the massive configuration. 

We now have a strategy for the constructing the particle phase space of a parton 
shower that involves massive particles. Starting with a system of 2 massless particles, 
construct a shower of massless vectors according to the procedure described above. 
In each antenna, let the momentum fractions of the (massless) final particles a, b be 
Wa, Wh- Rescale within the antenna by A and use the massless vectors and this value 
of A to compute the splitting probabilities. For example, for the splitting described 
in the previous section with particle a massive, the equation for A is 

^a + A(|4| + |4|) = Q. (34) 
The splitting probability is given by 



dProb = A^cT^(-77) / dwadwb ■ y ■ {^)r. ..^ , ^ tt r S . (35) 
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To evaluate the splitting function <S we need the flatted and sharped vectors a' and 
a^. The first of these is given by 

a^ = ^(£;„/A|4| + l)A4, (36) 

and a" is the massless vector of the same length pointing in the opposite direction. 
Once the configuration is chosen, the three new massless vectors are boosted back to 
the frame of the shower, and we are ready to generate the next antenna. When the 
shower is completed, the entire backbone must be rescaled to put the final massive 
particles on shell. In this prescription, the recoil due to emissions is done locally 
in each antenna to the extent that the particles are relativistic, but the recoil for 
nonrelativistic massive particles is distributed over the whole shower. 

There is one more complication that should be discussed. For a massless particle, 
the spin state is determined by the helicity in a way that is independent of frame. For 
a massive particle, a change of frame can rotate the spin. The helicity is preserved by 
rotations and by boosts along the direction of motion. Other boosts, at an angle to 
the direction of motion, change the spin orientation. In the massive particle shower 
described here, we ignore this effect. In any event, it is unimportant when the mas- 
sive particles are relativistic, and this accounts for most of the radiation from these 
particles. 



4 Antennae with a massive spin particle 

We are now ready to put together a catalogue of the antenna splitting functions 
that describe the emission of quarks and gluons in the showering of massive particles. 
We begin with the case of a spin massive particle S recoiling against a quark or a 
gluon. 

In the quark case, the antenna is described by an operator 

O^S'^{2\qR (37) 

where (2| is a spin-^ spurion that controls the quark polarization. Here and in the 
rest of the paper, we will analyze a subset of the various discrete choices from which 
the rest can be derived using the P and C symmetries of QCD. Here, for example, 
the two cases 

SQl SgLQL and Sq^ Sgnq^ (38) 

considered below suffice to provide all of the possble spin-dependent splitting functions 
for Sq Sgq and Sq Sgq. 
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The 2-particle matrix element of the operator (37) is 



M{0 ^ Sq^) = {2B) . 



(39) 



Then, for the 2-particle antenna Sq^ with S moving the 3 direction, 2 should be a 
massless fermion moving parallel to S. In the following, we will set 2 = A^. This choice 
follows the methods used in |H]. In that paper, the polarization vectors associated 
with operators O with nonzero spin are built from massless vectors 1 and 2, chosen 
in the directions of B and A, respectively. With this choice, the denominator of the 
expression ([T]) for the splitting function is again evaluated as (17). 



The 3-particle matrix elements of (37) are 

MiO^Sgnu) 



{A\h + c)ac) 
[cac) \bc\ 
{A^h){hac] 



[cac) (be) 



(40) 



Here again, we strip off the factors of g and color matrices. The final results are 
surprisingly compact. 

For an antenna containing a massive scalar and gluon, we need to find an operator 
that defines an antenna whose initial state includes a gluon of a definite polarization. 
For the antenna with a left-handed gluon, we may choose [8] 



= ^S^ (2|a-F|2) 



V2 



where 



a-F 



1 



-a cr 



(41) 



(42) 



This operator projects onto anti-self-dual gauge fields or left-handed physical gluons. 
The corresponding operator a ■ F can be used to define the antenna with an initial 
right-handed gluon. The two-particle matrix elements of (41) are 

M{O^SgL) = {2Bf M{O^SgR) = 0. (43) 

The zero for a gR is just as one should have expected. As above, we set 2 = Al . 



There are two types of 3-particle matrix elements of (41). First, the antenna can 
radiate a gluon. The corresponding matrix elements are 



MiO^SgLgL) 



1 

M 



[cac) 



[hah) 



12 



[cac) {be} 



KAU^ c ^ {A'c)\cah] 
M{0 ^ Sgngn) = , (44) 



following the pattern established in (40). Second, the gluon may split into a quark- 
antiquark pair. For this, we need the matrix elements 

{A^bf 



MiO^Sq^qL) 



{be) 



M{0^ Suqn) = . (45) 

The splitting functions derived from these matrix elements using ([T]) are listed sys- 
tematically in Appendix A. 



5 Antennae with a massive spin | particle 

In the same way, we can construct operators that correspond to the initial states 
of antennae involving a massive Dirac fermion Q with a quark or gluon. The massive 
fermion can have helicity ±|. Because the Q is massive, an initial left-handed Q can 
flip over after radiation to a right-handed Q, or vice versa. We have seen this already 
in the special case considered in Section 2. In this section, we will recall the results 
from Section 2 and compare them to those of the other three possible antennae of 
this type. 

The antennae with an initial state containing F and a quark can be arranged in 
a state with total spin about the axis of motion \J^\ equal to or 1. The spin case 
was considered in Section 2. The appropriate operator O is 

= QqR. (46) 

The matrix elements of this operator between two-particle Fq states are 

M{0^ QmL) = {A^B) M{0^ QL-qR) = (47) 

in our convention that should be used as the reference vector for Q. The three- 
particle matrix elements are then readily computed. If we use a" from the beginnning 
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as the reference vector for Q, (18) gives 

MiQigigL) = 



{ca^)[a^Qb) 
[a^c][cac) 
{a^b)[cQb) 
(be) [cac) 

{cJ)[a^Qb) 



'a^c]{a^a^)[cac) 



The antenna sphtting function can be constructed from these elements in the manner 
described in Section 2. 

The spin 1 case can be treated in the same way. As described in [8J and at the 
beginning of Section 4, we introduce hghthke vectors 1 and 2 in the direction of B 
and A, respectively. Then an appropriate operator to define this antenna is 

= Ql)[2qL. (49) 

The two-particle matrix elements of this operator are 

M{0 ^ QLqjt) = {Ah){2B) , M{0 ^ Qm^) = , (50) 

Thus, this operator does correctly represent the initial situation. We will set 2 = 
and 1 = i? in the following expressions. 

The splitting function for the antenna to radiate a gluon is computed from the 
three-particle matrix elements of this operator to Fgq final states. These are 



MiQLgRqR) 



{a^B)[A'{b + c)ac] 



(cac] (6c) 

,.,„ _ , m{a>B)\A*(h + c)ac\ 
= (a<a>>(c.c](fe) 



The splitting functions derived from these formulae and those in (48) are catalogued 
in Appendix A. 
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Figure 1: Feynman diagrams for the computation of the Fg — t- Fgg sphtting functions [8]. 



For the antennae with Q and a gluon, we again use the operator a ■ F to define 
the initial state as containing a gluon of definite left-handed polarization. There are 
two cases, with total spin | and |. For the spin ^ case, the appropriate operator is 

= -^Qa-F\2). (52) 

The dominant two-particle matrix element of this operator is 

MiO^QL9L) = {A'B){2B) . (53) 

If we recall that the vector 2 is identified with A^, we see that this puts the initial 
Q and g into just the correct orientation. The matrix elements to Qidn, QrQl^ and 
QrQr all vanish if 1 is taken parallel to B. 

The splitting functions for the radiation of a gluon from this antenna are given 



by the matrix elements of (53) to Qgg final states. As in pJJ and in [8], these matrix 
elements are given by the computation of the set of diagrams shown in Fig. [Tj The 
last diagram in the figure comes from the two-gluon vertex of the operator a ■ F. The 
third diagram is required to make the computation gauge- invariant. Its origin is most 
easily seen by thinking of the Q as a color octet. Then this diagram is obviously an 
essential contribution to the radiation from the Qg dipole. 



With this observation, we find for the three particle matrix elements of (53) 



[oc\[a*a"\ I [cac] {oab\ J 



{a'b){A'b){bac] 

MiQigRgL) = — — tjt-. — 
{cac\ {be) 

{a'c){A'c){cab] 
M{Ql9r9r) = 
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MiQugngL) 
MiQngLgR) 



m 



m{a%){A%){hac] 



(atta^) (cac] (6c) 
m(a«c)(A^c)(ca6] 

-M(QR5R5i?)=0. (54) 



The case of a antennae in the spin | state is treated similarly. The operator 
that defines the initial state is 

= -^Ql]{2\a.F\2) . (55) 

The two-particle matrix elements of this operator are 

MiO^QRgL) = {lA'){2Bf (56) 

and all other matrix elements are equal to zero for the choice of 1 parallel to B. We 
will set 2 = and 1 = i? in the expressions that follow. 

The three-particle matrix elements of ( [55| to Qgg final states are 

_ [a'B] f {A'b){A\b + c)ac) {A'c){A\b + c)ab) 

MiQngLg^) = + 



fa'') [6c] I (cac] {bab] 



m{A%f {[a^B]{bac] + [ah]{bcB]) 
[afa^] (cac] (6c) 

m{A^cf ([a«B](ca6] + [a»6](c65]) 



MiQugR-gR) = 

, mlJB] ((A^b)(A\b + c]ac) (#c) (#(6 + c)a6) 

MiQigRgL) = 

MiQLgigR) = r . w//, \ 

[afa"] (6a6] (6c) 

A^(QL^?i?^i?) = 0. (57) 

The splitting functions for Qg — )■ Qgg that are derived from these expressions and 
those in (54) are catalogued in Appendix A. 
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The Qg antennae can also radiate by gluon splitting to a pair of quarks. For the 
spin I case, the relevant matrix elements are 

MiQLqRQL) = 

{a^c){A'c) 
{be) 
m{a%){A^h) 
{aia}>){hc) 
m{Jc){A^c) 
(atta^)(6c) 

For the spin | case, the matrix elements are 

^^.^ , [a'B]{A%f 
MiQuq^qL) = (^cT^ 

MiQuq^qR) = - ' 



MiQLqLQR) = 
MiQnqRqL) = 
MiQmLqR) = 



MiQLqRqh) 



{be) 

m[a^B]{A^b) 
[Ja^]{bc) 



(58) 



x.fn - ^ m[JB]{A'c)^ , , 

The splitting functions for Qg — )■ Qgg that are derived from these expressions are 
catalogued in Appendix A. 



6 Antennae of a pair of massive particles 

After a pair of massive scalars or fermions are produced, their first emission of 
a gluon is described by an antenna in which the two massive particles both appear. 
For a complete description, we need the splitting functions for these antenna as well. 
These formulae are somewhat more comphcated than those derived above, since some 
of the simplifications that are possible when the particle b is massless no longer apply. 
There is little additional complexity in the cases in which the two massive particles 
have different masses, so we will write the formulae for that more general situation. 



17 



The case of a pair of scalars is relatively straightforward. The scalar particles 
themselves are spinless, so there is only one case, described by the spin-0 operator 

O = SlS2 . (60) 

The matrix element of this operator to create the state S1S2 is simply 1. The matrix 
elements for gluon emission are 



[a"c] I [cac) {cbc] 

KA(Q Q\ ^ /^^^^^] <^'^^U ^a^^ 

-A^(5i^?^5,) = -^|^-^|. (61) 

Each expression can be brought down to one term using the Schouten identity 

(ca/] {dhg] - {daf] {cbg] = - {cd) [fahg] . (62) 

This identity is valid when a and h are massive vectors, possibly with different masses; 
c, d, /, and g must be massless. To prove the identity, write a as a linear combination 
of and a". Using (62), 



(cac] [chc) 

M{S,g,S,) = . (63) 

[cac) {cbc\ 

The splitting functions are readily assembled from these expressions. 

For the antenna of a massive fermion and a massive scalar, the general case is 
described by the spin | operator 

= Q, l)S2 . (64) 

The two-body matrix elements of this operator are 

A^(QiL^2) = {A'l) (65) 

and zero for Qir. If we take 1 = following the prescriptions above, 

\{A'l)\' = {E, + K){E2 + K) , (66) 

where i?2, and K are the two energies and the momentum in the antenna center 
of mass frame. 
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The matrix elements for the operator (64) to create QgS states is given by the 
expression 



M 



V2 



u[a) 



+ 4 + m) 



2b ■ e(c) 



(67) 



[cac) ' ' \cbc) 

where e(c) is the polarization vector of the gluon. A convenient way to treat this is 
to manipulate 

/(c)(^ + ^ + m) = 2a-e(c) + /(c)^ (68) 
that gives zero when applied to u[a). The first 



plus a term proportional to 



term in (68) combines with the last term in (67) to give an amplitude proportional 



that of the scalar-scalar case, (61) or (63) above. The term with ^(c) vanishes for gji 
and gives a simple but nonzero term for gi^. The final results for the two amplitudes, 
after dropping the factor of {gT°'), are 



MiQLgRS) 
MiQRgRS) 
MiQLgLS) 
MiQRgLS) 



{ah){cabc) 
{cac] [cbc) 
mi{a^l) (cabc) 
{a^a^){cac\[cbc) 
{a^l)[cabc] (a^c)(cl) 



+ 



[cac) {cbc\ [cac) 



mi 



{ah)[cabc] _ (a«c)(cl) 
[cac) {cbc\ [cac) 



(69) 



Here mi is the mass of the fermion Qi. The formulae apply for any values of the 
masses of the fermion and scalar, as long as the 4-vectors a and b are properly on 
mass shell. 



The decomposition of the gluon coupling to a massive fermion given in (68) is 



equivalent the representation of this coupling by the second-order Dirac equation, in 
which the fermion is replaced by a field with a scalar-type coupling and a magnetic 
moment coupling. The single-gluon magnetic moment coupling has a chiral structure 
and vanishes for specific combinations of the fermion and gluon spin. This second- 
order Dirac formalism is discussed in more detail in 



For massive fermions, there are two cases, corresponding to total spin and 1 
along the antenna axis. For the spin case, we could use the operator QlQl to 
create the antenna, similarly to the choices in Sections 2 and 5. However, in the case 
in which both fermions are massive, that operator creates both QlQl ^^cl QrQr 
states. We will avoid that problem here by taking the operator that creates an initial 
state of QlQl to be 

= Qil){2Q2 (70) 
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Figure 2: The single Feynman diagram for the computation of the gg — )• gFF spHtting 
function [Hj. 



The two-body matrix elements of this operator are 



M{Q^lQ2l) = {A'1){2B') 



(71) 



and zero for the other three hehcity states. Similarly, for the spin 1 case, we will use 
the operator 



to create an initial state of QlQr- The two-body matrix elements of this operator 
are 



and zero for the other three helicity states. The QgQ matrix elements of these 
operators are easily computed using the methods presented earlier in this section. 
The results for the splitting functions are tabulated in Appendix A. 



7 Antennae with massive particle production 

There is one more situation that we must consider. At very high energies, massive 
particles can be produced by gluon splitting. At the LHC, for example, parton-parton 
scattering can give quark-gluon and gluon-gluon collisions with center of mass energies 
well above 1 TeV. Final state gluon antennae in these coUsions can produce pairs of 
top quarks. The pair production amplitudes are relatively simple, since each requires 
only one Feynman diagram, as shown in Fig. [2] for the gg — )■ gtt case. The final pair 
of heavy particles must have equal mass and equal spin. However, there are a large 
number of cases to enumerate. The massive scalar or fermion pair can be formed 
from a spin ^ or a spin | qg antenna or from a spin or spin 2 gg antenna. 

For scalar pair production, the formalism is actually quite simple. The spin | and 
spin I qg antennae can be represented by the operators 



= Q^l)[2Q2 . 



(72) 




(73) 
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03/2 = -^gl](2|a-F|2) . (74) 



If the gluon splits to a pair of scalars, both cases involve the operator a ■ F dotted 
with the gSS vertex. This product is 

\[{H4)l'-l'{Hmh-c),^[h,c] (75) 
so that the relevant three-particle matrix elements are 

M{<i,-SS) = -^^^ (76) 

Sbc 



for the spin | case and 



for the spin | case. 



She 



Similarly, the spin and spin 2 gg antennae, corresponding to the gLdL and gngL 
initial states, can be represented by the operators 

Oo = ltT[{w . FY] 

02 = [l\a-F\l]{2\a-F\2) . (78) 



The manipulation (75) again gives a simple form for the three-particle matrix ele- 
ments. The corresponding splitting functions are given in Appendix A. 

For the case of massive fermion pair production, this formalism is necessarily more 
complex. With the choice of helicity states that we have used throughout this paper, 
the vertex to create a pair of massive fermions is a Dirac matrix. For the case of a 
final-state QrQl, for example, this matrix has the form 

Then the matrix element of a ■ F contains the structure 

{Ri[{b + c),V]R2) (80) 

with a commutator bracketed between reference vectors -Ri and R2. However, the 
frame-dependent choice of the vectors b^, c" makes it difficult to simplify this expres- 
sion further. It is true that {b+c) = Q — a, where a is now massless. In some cases, we 
have Ri = a, in which case the a term cancels. In other cases, we have Ri = 2 = A\ 
so that the a term vanishes if a is coUinear with A. We list the full expressions for 
these splitting functions in Appendix A. 
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8 Conclusion 



In this paper, we have provided new materials for the construction of parton 
showers that include massive spin and spin | particles. Wc hope that this formalism 
we have presented will be useful in describing the QCD dynamics of the top quark 
and other heavy particles at LHC. 



A Catalogue of massive antenna splitting functions 

In this appendix, we catalogue the various antenna splitting functions for massless 
particles derived in this paper. Antenna splitting functions not listed are equal to 
cases listed below that are related by the P and C symmetries of QCD. 

A.l Splitting functions with one massive scalar 

A. 1.1 Spin I antenna: initial dipole Sq^ 



Q_ 
2K 



{A\b + c)ac) 



[cac) [be] 



Q {A^b){bae] 



2K [eae){he) 



(81) 



A. 1.2 Spin 1 antenna: initial dipole SgL 



S{SgLgL) 
Si^SgngL) 



{2Ky 



[6c] L [cac) 



[hah) 



S{SgLgR) 
SiSgRgn) = 

^{SqjiqL) 



1 


{A^h)\hae] 2 


(2X)2 


[eae) (he) 


1 


{A^e)\eah] 2 


(2X)2 


[hah) (he) 







1 


[bA^b)^ 


(2i^)2 


[heh) 
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S{SqLqR) 



{2KY [hch) 



(82) 



Note that the last two expressions are aheady squared and evaluate to values that 
are real and positive for the case of a final-state antenna. For example, 



1 {2b ■ A^] 
{2Kf 2b ■ c 



A. 2 Splitting functions with one massive fermion 
A.2.1 Spin antenna: initial dipole Ql^l 



SiQLQLqL) 

SiQLgmL) 
S{Qr9rQl) 



Q_ 

2K 

Q_ 
2K 
rri^Q 



[a^c] [cac) 
{a^b)[cQb) 



{be) [cac) 
(a^c) [a- 



2K 



{a^a^)[a^c][ 



cac) 





{a%)[cQh) 


2K 


{a^a^){hc) [cac) 



(83) 



A.2.2 Spin I antenna: initial dipole QlQl 



S{Ql9l9l) = 

SiQLgRQL) = 

S{Ql9l9r) = 
S{Ql9r9r) = 
S{Qr9l9l) = 



1 


{A^h){Q^[a^ac) - 


(2ir)2 


[a^a^] {cac 


1 


{a%){A^b){bac] 2 


(2ir)2 


(cac] {be) 


1 


{a^c){A^c){cab] ^ 


{2KY 


{bab]{bc) 



+ 



[aia^]{bab][bc] 



m 



{2K\ 



{A%){{ahQc) - Q'(a«c)) {A^c){{ahQb) - Q^{a^b)) 



(atla'')(cac][6c] 



+ 



{a^a^){bab][bc] 
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SiQRgLgR) 
SiQugRgR) 
S{QLqRqL) 

SiQRqRqL) 
SiQRqLqn) 



m 



{2KY 

2 

m 



{a%){A%){hac] 



{a^a^){cac]{hc) 
{ah){A^c){cah\ 



{a^a^){hah]{hc) 



{2Kf 


1 [ha}'h)[hA^h) 
(2KY \bcb) 
1 [ca!'c)[cA^c) 

rri^ [ha%)[hA^h) 
(2KY [aiaai)[hch) 
[ca^c)[cA^c) 
(2K)2 [aiaai)[bcb) 



(84) 



As in A. 1.2, the last four expressions here are already squared and evaluate to 
real, positive values. 



A. 2. 3 Spin 1 antenna: initial dipole QhqR 



SiQLgRqR) 

SiQigLqR) 
S{QRgRqR) 
SiQRgmR) 



1 


{a^B)[A\h + c)ac] 


(2^)2 


(cac] {he) 



{2Ky 



[A^b]{[aiac)[bQB) +m\cB)[a% 



[ah^]{cac][bc] 



2 

m 


{a^B)[A'{b + c)ac] 


2 


{2KY 


{a^a}'){cac]{bc) 




2 

m 


[A^b]{{a^B){cab] + 


{ah) {Bob]) 2 


{2KY 


{a^o!'){cac]{bc) 



(85) 



A.2.4 Spin | antenna: initial dipole QRgL 



SiQLgLgh) = 





[a^B] 


(2QX)3 


[a^a)>][bc] 



(cac] 



+ 



{bab] 
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S{Ql9r9l) 

S{Ql9l9r) 
S{Ql9r9r) 
S{Qr9l9l) 

S{Qr9r9l) 

S{Qr9l9r) 
S{Qr9r9r) 

^iQLQRQL) 
SiQLQLQR) 

SiQRQRqL) 
S{QRqL<lR) 



{2QKY 



{A%) {[a^B]{hac] + [ah]{hcB]) 



{2QKY 




[a^a^]{cac]{hc) 
{A^cf{[a^B]{cah] + [a%{chB]) 







{2QKY 


[he] 



[a^a!>]{hah]{hc) 



(cac] 



+ 



[hah] 



{2QKY 



1%) 



[a'c]{hQB]+m'^^[cB] 



{2QKf 



(cac] (6c) 



{bah]{bc) 



(a«a^) 



= 



m2Q2 


[5a«5)[6#6)2 


(2Q;^)3 


[a^aa^)[hch) 




[Ba^B)[cA^cf 


{2QKf 


[flttaatt) [6c6) 




[Ba''B)[bA%f 


{2QKf 


[bob) 




[Ba^B)[cA^c)'^ 


{2QKY 


[bob) 



(86) 



As in A. 1.2, the last four expressions here are aheady squared and evaluate to 
real, positive values. 



A. 3 Splitting functions with two massive scalars 
A.3.1 Spin antenna: initial dipole S1S2 



(cabc) 



{cac] [cbc) 
[cabc] 



[cac) {cbc] 



(87) 
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A.4 Splitting functions with a massive fermion and a massive scaleir 

A.4.1 Spin I antenna: initial dipole Q1S2 



SiQlL9LS2) 
S{QlRgLS2) 

S{Qil9rS2) 
S{QiLgRS2) 





{a^B^){cabc) 


{Ei + K){E2 + K) 


[cac] [cbc) 



mlQ^ 


1 


{Ei + K){E2 + K) 





(cac] 

{a^B^){cabc) _ {ah){cB')^ 2 
(cac] [chc) (cac] 





{a^B^)[cabc] 2 


{Ei + K){E2 + K) 


[cac) (cbc] 



{Ei + K){E2 + K) 



1 {a^B^)[cabc] 



(a'ia^) [cac) {cbc] 



(88) 



A. 5 Splitting functions with two massive fermions 

A. 5.1 Spin antenna: initial dipole Q\lQ2L 



S{Qil9lQ2l) = 



S{QiLgLQ2R) 



S{QiRgLQ2R) 





{a^B^){cabc){A^b^) 


{{E^ + K){E2 + K)y 


(cac] [cbc) 



(a^c) (cS^) {A%^) {a^B^) (A^c) {cb^) 



{(E, + K){E2 + K))^ 



{cac] 

1 ({a^B)^{cabc){A^b^) 



[cbc) 



{b^bi){ {cac][cbc) 

{a^c){cB^){A^b^) {a^B^){A^c){cb^)^ 2 



S{QiRgLQ2L) = 



m2Q2 



{{E, + K){E2 + K)y 



{cac\ 

1 f{a^B)\cabc){A^b^) 



[cbc) 



(a^a^) I {cac][cbc) 

{a^c){cB^){A^b^) {a^B^){A\:){cV)^ 2 





1 


{{E^ + K){E2 + KW 


(atta^)(6fl6^) 



(cac] [c6c) 
(a«5^)(ca6c)(#6«) 



(cac] [c6c) 

(a«c)(c5^)(#6«) (a«5^)(#c)(c6«)^ 2 



(cac] 



[c6c) 
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S{QiLgRQ2L) 
S{QiLgRQ2R) 
S{QiRgRQ2L) 
S{QiRgRQ2R) 





{a}'B^)[cabc]{A^y') 


{{Ei + K){E2 + K))^ 


(cac] [cbc) 



{{Ei + K){E2 + K)y 
{{Ei + K){E2 + K))^ 
{{Ei + K){E2 + K)y 



1 

1 {a^B^)[cabc]{A^bi) 



i&bi) {cac][cbc) 

1 {a^B)^[cabc]{A^b^) 
(a^a^) (cac] [cbc) 

1 (a»B^)[ca6c](#6«) ^ 
JaMyim^ {cac][cbc) 



(89) 



A. 5. 2 Spin 1 antenna: initial dipole Q1LQ2R 



S{QiLgLQ2L 

S{QlLgLQ2R 
S{QlRgLQ2L 

S{QiRgLQ2R 

S{Q\LgRQ2L 

S{QiLgRQ2R 
S{QiRgRQ2L 
S{QiRgRQ2R 



{El + K){E2 + K))^ 



1 r ia^B^) (cabc) [A^b^] {a^c) {cB^)[A^b^] \ ^ 



(El + K){E2 + K))^ 



(6^6«)l {cac][cbc) {cac] j 

(a^B^) (cabc) [A^b^] {a^c) (cB^) [A%^] \ 2 



{El + K){E2 + K))^ 



cac] [cbc) [cac 

1 r ( a^B^) (cabc) [A^^] {a^c) (cB^) [A^^] j ^ 



{Ei + K){E2 + K))^ 



(atta^) I {cac][cbc) {cac^ 
1 ((a^B^){cabc)[A^b^] {ak){cB^){A^b')^ 2 



{Ei + K){E2 + K)y 



(a^a^) I {cac\[cbc) {cac^ 
1 ({a^B^)[cabc][A%i] {a^ B^)[A^ c][cb^^ ^ ^ 



{Ei + K){E2 + K))^ 
m\m\Q'^ 



[6^6tt] I {cac][cbc) [cbc) 

{(eB^)[('(iiH]lV}y] {(e B')[A\][cb^] 2 



{Ei + K){E2 + K)Y 



(cac] [cbc) [c6c) 

1 r (a«S^) [cafec] [A^6^] (a^S^) [A^c] [c&«] | 



(atta^)[6^6»] 



(cac] [c6c) 



[c6c) 



(£;i + x)(£;2 + i^))^ 



/ (a«E^)[ca6c][A^&^] _ (a^E^)[#c][cy] ) 



(a*a'') I (cac] [c6c) 



[c6c) 



(90) 



A. 6 Splitting functions with pair production of scalairs 

A.Q.I Spin antenna: initial dipole gigL 



S{gLSS) = — 



{a[b. c]a) 



Sbc 



(91) 
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A.6.2 Spin \ antenna: initial dipole qidL 



{a[b, c]A) 



Sbc 



A. 6. 3 Spin | antenna: initial dipole gnQL 



SiqnSS) = 



\aB\^A\b,c\A) 



Sbc 



A.6.4 Spin 2 antenna: initial dipole gnQL 



S{gRSS) = ^ 



[aBf{A[b,c]A) 



Sbc 



A. 7 Splitting functions with pair production of fermions 

A. 7.1 Spin antenna: initial dipole QlQl 



S{gLQLQL) 
SigLQLQn) 
SigLQnQL) 
S{9lQrQr) 



2 

m 


{aQb^]{c^a) 


Q'4c 


[b^bi] 



Q's' 



be 



{aQc^]{b^a) + 



2 



he 



[6^6«](c^ctt) 



m 



(6^6«)[c^cfl] 



(aQ6«](c«a) 



Q's' 



be 



{aQc^]{b^a) ^ {aQb^]{c^a) ^ 



[b^bi] 



A. 7. 2 Spin \ antenna: initial dipole QlOl 



SiQLQLQL) 



2 

m 


(aQb^ic^A) 




[b^b^ 



+ 
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^ {A{Q-ap]{&a) ^ {A{Q - a)c«](6^a\ 2 



[c^c«] 



be 



{aQc^]{b^A) + 



m 



{aQh^]{c^A) 



+ {A{Q-a)c^]{}}'a) + 



m 



{A{Q-a)h^]{c^a) 



SiqLQRQh) 



be 



aQ})']{c^A) + 



m 



{aQc^\{h^A) 



+{A{Q - a)})>]{c'a) + j^^^^^{A{Q - a)c«](6«a) 

{aQh^]{c^A) 

(A(Q-a)c^](6"a) {A{Q - a)l)']{c^a': ^ 
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{aQc^]{h^A) 
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A. 7.3 SpiiJ I antenna: initial dipole qngi 



SiqRQiQL) 
SiqnQLQR) 
SiquQuQi) 
SiquQuQu) 



m 



be 



[aB] 



{A{Q-a)b^{c^A) ^ {A(Q-a)cmb'A) \ 



be 



be 



[b^bi] 

[aB]l^{A{Q-a)c^]{b^A) + 
[aB]\{A{Q-a)b%^A) + 



7 



m 



[6^6tt](c^c«) 



m 



{b^U)[&d] 



(A(Q-a)fc«](c«A)| 



{A{Q - a)c' 



m 



be 



[aB][ 



{A{Q - a)c^]{bU) {A{Q-a)b'']{cU) \ ^ 
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A. 7.4 Spin 2 antenna: initial dipole gnQL 



SigRQLQL) 
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Q'4e 
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[aB] 



. r (A(g - a)6«] (c^A) (A(g - a)c«] (6^ A) 1 



[b^bi] 



+ 



be 



aBY{{A{Q-a)c^]{b^A) 



m 



[b^U]{&c^) 



{A{Q-ap]{c^A) 
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SigRQnQL) 



he 



[aBf\{A{Q-a)})']{c^A) + 



m 



(A(Q-a)c«](6»A) 



m 



be 



{A{Q-a)c']{biA) {A{Q-a)b']{c^A) \ 



(9^ 



B Spin-dependent Altarelli-Parisi functions for massive par- 
ticles 

In this Appendix, we present the spin-dependent Altarelh-Parisi sphtting functions 
for massless and massive particles. The massless cases were derived in the original 
paper of Altarelli and Parisi [3J . Spin-summed Altarelli-Parisi functions for the cases 
with massive particles arise in NLO QCD calculations for supersymmetric particle 
production. They have been catalogued by Catani, Dittmaier, and Trocsanyi in [20] . 
The spin-dependent functions can be worked out by textbook methods. Here we 
present these functions in a representation convenient for comparison to the antenna 
splitting functions derived in this paper. We omit the overall color factor of Nc and 
divide by 2 so that the splitting accounts the contents of an individual antenna. 

Note that, since we work at the leading order in Nc and normalize to a single 
antenna, there is no difference between the splitting function for a heavy quark or a 
gluino to radiate a gluon. Thus, there are only two cases, the cases of a heavy scalar 
S* or a heavy quark Q radiating a gluon. The cases of a heavy particle splitting to a 
heavy particle by radiating a gluon are given by the same expressions with z — )• {1—z). 

For S — )■ gS, 

PiS ^ SgLS) 
P{S ^ SgnS) 



1 — z 

Pj. + z'^Trfi z 

pip 1 — z 
Pj, + z'^m? z 



(99) 



For Q gQ, 

P{Ql ^ QigL) 



Pj, + z^m^ z 

Pt (1 - 
Pp + z^m^ z 
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.2 „4 



P{Ql Qr9l) = 2 ^2 2 ~ 

P{Ql ^ Qr9r) = (100) 
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